Optical ratchets with discrete cavity sohtons 
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We propose a setup to observe soliton ratchet effects using discrete cavity solitons in a one-dimensional array 
of coupled waveguide optical resonators. The net motion of solitons can be generated by an adiabatic shaking 
of the holding beam with zero average inclination angle. The resulting soliton velocity can be controlled by 
different parameters of the holding beam. © 2008 Optical Society of America 
OCIS codes: 190.0190, 190.5530, 230.1150 



Spatial solitons in optical nonlinear cavities have at- 
tracted a lot of attention during last decades, see, e.g., 
reviews in Refs. 1,2 and references therein. They rep- 
resent nondispersive localized structures supported by 
a fiat holding beam (pump) and result from the inter- 
play between nonlinearity and diffraction, as well as gain, 
losses and internal feedback of the system. Due to mul- 
tiple reflections of light at the boundary mirrors, which 
form the cavity, light interaction with the nonlinear ma- 
terial inside the cavity is effectively increased. As a con- 
sequence, spatial cavity solitons (CSs) may be formed at 
essentially reduced input powers, as compared to con- 
ventional spatial solitons in waveguides. All the param- 
eters of cavity solitons (including energy and phase) are 
completely determined by the parameters of the hold- 
ing beam, which adds more flexibility to the control of 
the process of creation and annihilation of CSs as well 
as to their evolution. Together with a high robustness, 
usually inherent for solitons, the above properties turn 
CSs into good candidates for information storage units 
in all-optical devices. Cavity solitons have been experi- 
mentally observed in semiconductor microcavitics^ and 
single- mirror feedback loops with nonlinear elements.'^ 

The presence of gain and losses in the system intro- 
duces several qualitatively new features to the soliton 
dynamics. In particular, a unidirectional soliton motion 
under the influence of AC zero-mean external forces can 
be realized, provided certain symmetries of the system 
are broken.^' ^ Such a soliton ratchet effect has been 
successfully observed in experiments with an annular 
Josephson junction, revealing a unidirectional topolog- 
ical fluxon motion driven by biharmonic microwaves of 
zero average."^ Up to now the soliton ratchet effect has 
been discussed only for solitons bearing a non-zero topo- 
logical charge. In such a case the excitation can not be 
destroyed due to large energy barriers in the system. 
In this Letter we extend the concept of soliton ratch- 
ets to much more fragile optical cavity solitons having 
zero topological charge. We consider a nonlinear optical 
waveguide array with dielectric mirrors at the end facets 
(array of coupled zero-dimensional optical resonators), 
driven by an external field (pump), see Fig.^ Recently 
it was demonstrated that discrete cavity solitons (DCSs) 
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Fig. 1. Schematic setup of coupled waveguides with mir- 
rors at the end facets. 



can be excited in such a system within a reasonably wide 
range of control parameters.® The evolution of the trans- 
mitted field pattern at the output facet obeys a discrete 
nonlinear Schrocdinger (DNLS) type equation:® 
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where A„ and A™ correspond to the renormalized ampli- 
tudes of transmitted and input fields in the nth waveg- 
uide, respectively. The time r t/To is measured in 
terms of the photon lifetime inside the cavity Tq = 
2d/[vg{l — p^)], which may be of the order of picosec- 
onds or larger. Here d is the length of each waveguide, 
Vg is the group velocity of light inside the waveguides and 
p is the real-valued reflection coefficient of the mirrors, 
(1 — p^) ^ 1. The effective damping coefficient is rescaled 
to unity in Eq. while parameters C = Co/(l — p^) 
and A = Ao/(l — p^) + 2C define the effective cou- 
pling between adjacent waveguides and detuning from 
the linear resonance, respectively (Aq = 2/3od, where 
is the wave number corresponding to the laser fre- 
quency uJo). Both parameters can vary in a wide range 
by adjusting the reflectivity of mirrors p, as well as the 
frequency of incident light and the distance between ad- 
jacent waveguides (the latter determines the value of 
Co). The nonlinear Kerr coefficient a can be rescaled 
to a = ±1. Increasing the coupling parameter C the 
model in Eq. asymptotically approaches the continu- 
ous Lugiato-Levefer model^ for a one-dimensional res- 
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Fig. 2. Dark (a) and bright (b) DCSs profiles. Black lines 
and circles: stationary DCSs for the input power |ap = 
3.3 at normal incidence 4>in = 0. Model parameters are: 
(a) a = -1, A = 3, C = 0.55; (b) a = 1, A = -3, C = 
15. Gray lines and crosses: snapshots of moving DCSs for 
the ratchet case with the same parameters as in Fig. 

onator based on a homogeneous-core slab waveguide. 
Then the discrete Laplacian in Eq. ^ is replaced by 
the second spatial derivative d'^/dx^ of a continuous field 
A{x,t). 

For a constant input field 
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4>in determines the incident angle of the pump light. Nor- 
mal incidence (j)^ = Eq. supports stationary DCS 
solutions^ in a certain regime of the input power |ap and 
detuning A, sec Fig.[5| A tilted holding beam with the 
incidence angle 4>in ^ induces a transverse force acting 
on the soliton. It results in a moving soliton which can 
be observed both in continuous^*^ and discrctc^^ models. 

Let us now consider the case, when the incidence angle 
of the holding beam is adiabatically slowly oscillating in 
time with zero mean: 0m (r + T) = (j)in{T) (T ^ 1) , 
/J" 4>in{T)dT ~ 0. We want to study whether the cor- 
responding induced AC zero mean force can result in a 
unidirectional (in average) ratchet motion of the soliton, 
disregardless the initial conditions (i.e. the initial phase 
of <f>in{T), position of the soliton, etc.). For the period of 
the holding beam oscillations T being much larger than 
the characteristic photon lifetime in the cavity, T ^ 1, 
the soliton dynamics is locked to the beam phase, i.e. 
there exist only a unique attractor solution. As a first 
step we perform a symmetry analysis^ of Eq. Sup- 
pose we identify a symmetry operation 5, which involves 
shifting of spatial and temporal coordinates and/or mir- 



ror reflections, leaves the dynamical equations invariant 
and, at the same time, changes the sign of the soliton ve- 
locity. If applied to any trajectory, it will again generate 
a trajectory. Since the attractor is unique, the symmetry 
operation applied to the attractor will map it onto itself, 
and the average velocity on this attractor will be exactly 
zero, so that the soliton will perform periodic oscillations 
in space and time. Therefore, breaking such symmetries 
of the underlying model equations is a necessary condi- 
tion to observe a soliton ratchet effect. 

Far away from the soliton center the field distribution 
asymptotically approaches the homogeneous (in inten- 
sity) ground state. The corresponding ground state field 
distribution A^^{t) w a'^'^(T) exp(i(/)*"n) and its inten- 
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P are obtained from^ 
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Then, we can define the position of the center of soliton 
X{t) and its velocity V{t) 

""^"^ = EJKMP-/-(r)| ' = • 

(4) 

Now we note, that any single-harmonic variation of 
4'in{T) = sin[cij(T — To)] posscsscs thc time-shift sym- 
metry 

0„ (T + r/2) = -0,„(T) . (5) 

Provided this symmetry holds, the dynamical equation 
^ remains invariant under the operation Sa 
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At the same time, this operation changes the sign of the 
soliton velocity. Therefore, no rectification of a single- 
harmonic AC force is possible with cavity solitons. How- 
ever, the above symmetry can be broken e.g. by choosing 
a biharmonic variation of thc incidence angle: 



(/'jti(t) = 0a sin(wT) -f sin(2cjT + 6) 
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Yet another possibility is to use two holding beams, 
which oscillate at the first and second harmonic fre- 
quency, respectively: 

A- = a exp (^)«] + ^ exp [z^^^^ (t)^] , (8) 

"/"in = sin(wr)n , <j)f^ = 0b sin(2LJT + e)n . (9) 

In both cases a unidirectional propagation of DCSs can 
be observed, despite the fact that thc averaged value of 
the incidence angle (and thus that of the force acting 
on thc DCS) is zero, sec Fig. 13 Thc resulting average 
velocity of thc soliton strongly depends on the system 
parameters and parameters of shaking beam(s) [cf. max- 
imum velocities of dark and bright DCSs for different 
holding beams and coupling parameters in Fig. ^ . One 
of the possible ways to control the DCS motion is to 
adjust the relative phase 9 in 0,©, see Fig. ^ 
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Fig. 3. Density plots of |A„p for soliton motion in 
the ratchet case with initial conditions corresponding 
to stationary DCSs. Model parameters are the same as 
in Fig. 121 (a) Dark DCS driven by the single shaking 
beam ©, 0, |ap = 3.3, 0„ = 0,, = 0.025, u = 0.01, 
= 1.72 • TT. (b) Bright DCS driven by two shaking 
beams ®, |ap = \b\^ = 1.2, 0„ 0^ = 0.001, uj = 0.002, 
6' = 0.57-7r. 

For slowly shaking holding beams the ratchet effect 
can be estimated through the DCS velocities as functions 
of constant tilt V(</'m). The resulting average velocity of 
a DCS can be obtained as 

V = ^ £ V[c|^^n{r)]dT (10) 

We note, that in the continuum limit of Eq. ^ the soli- 
ton velocity V{(j)in) is a linear function of the incidence 
angle. Thus, for any choice of periodic functions (/)i„ 
with zero mean the resulting average velocity l|10|l is zero 
in the adiabatic limit. The ratchet effect for slowly vary- 
ing incidence angle(s) (pin can appear only for nonlin- 
ear functions V{(f>in)- This nonlinearity is induced by the 
discreteness of the system due to Peierls-Nabarro poten- 
tial.^^ In that case the resulting soliton net motion does 
not depend on the actual choice of the frequency lu which 
can be taken arbitrary small. 

To conclude, we have demonstrated the possibility to 
generate a rectified motion of discrete cavity solitons by 
means of shaking holding beams with zero average in- 
clination. The necessary condition for the observation of 
this soliton ratchet effect is the violation of the symmetry 
Sa © of the underlying model equations. This can be 
done e.g. by bi-harmonic variation of the holding beam 
inclination angle or by a superposition of two shaking 
holding beams with different frequencies. The velocity 
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Fig. 4. Soliton displacement per period T as a function 
of the relative phase 9 for dark (green line) and bright 
(black line) solitons. All the parameters are the same as 
in Fig. 121(a) and (b), respectively. 

of the resulting soliton net motion can be adjusted by 
the parameters of the holding beams (w, 0a, 0ti, ^)- This 
opens a new prominent perspective for soliton steering 
and various all-optical switching schemes. Our results 
are also instructive for the general problem of nontopo- 
logical soliton ratchets in spatially extended discrete and 
continuous systems. 
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